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Abstract 

^ I The key question in the interaction of antinucleons in the nuclear medium concerns the 
■ deepness of the antinucleon-nucleus optical potential. In this work we study this task in 
the framework of the non-linear derivative (NLD) model which describes consistently bulk 
properties of nuclear matter and Dirac phenomenology of nucleon-nucleus interactions. 
We apply the NLD model to antinucleon interactions in nuclear matter and find a strong 
decrease of the vector and scalar self-energies in energy and density and thus a strong 
^ suppression of the optical potential at zero momentum and, in particular, at FAIR energies. 

^ ■ This is in agreement with available empirical information and, therefore, resolves the issue 
^ concerning the incompatibility of G-parity arguments in relativistic mean-field (RMF) 

models. We conclude the relevance of our results for the future activities at FAIR. 

(N 

^ ■ Keywords: relativistic hadrodynamics, non-linear derivative model, nuclear matter, 

Schrodinger equivalent optical potential, proton-nucleus optical potential, 
Tj- ' antiproton-nucleus optical potential 

in 

O ■ 1. Introduction 

The in-medium nucleon-nucleon interaction has been an object of intensive theoretical 
and experimental research of modern nuclear physics over the last few decades, see for a 
^ [ review [if. The main finding was a softening of the nuclear equation of state at densities 
reached in intermediate energy nucleus-nucleus collisions, which was consistent with a 
variety of phenomenological |2| and microscopic [3] models. In addition the empirical 
saturation of the proton-nucleus optical potential turned out to be consistent with heavy- 
ion theoretical studies 0]. 

While the bare antinucleon-nucleon (NN) interaction has been actively studied, see 
Refs. ^ and references therein, empirical information on the in-medium interactions of 
antinucleons is still very poor. Antiproton production has been investigated theoretically 
in reactions induced by protons [gI and heavy ions in the SIS-energy region [3], where some 
data on antiprotons were available. Complementary studies of antiproton annihilation in 
nuclei [8|] and antiprotonic atoms [9] provided further insight on the optical potential at 
very low energies, however, with rather big uncertainties in the nuclear interior due to the 
strong annihilation cross section at the surface of the nucleus. 

In the near future the FAIR facility intends to study the still controversial and em- 
pirically less known high energy domain of the (anti)nuclear interactions in more details 
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than before. For instance, the nuclear equation of state for strangeness degrees of freedom 
and also the in-medium antinucleon-nucleon interaction are some of the key projects l3 • 



They are relevant for the formation of exotic (anti)matter systems such as double-strange 
hypernuclei and A-hypernuclei in antiproton-induced reactions in the PANDA experiment 
at FAIR pj]. _ 

The microscopic Brueckner-Hartree-Fock calculations of the in-medium A^A^-scattering 



have been carried out in [12| . On the other hand, a complementary theoretical background 
for phenomenological models builds the relativistic hadrodynamics (RHD). It is based on 
the relativistic mean-field (RMF) theory, which is a well established tool for infinite and 
finite nuclear systems 13 1. However, as already shown many years ago 0], there are still 
unresolved problems in RMF models, when applying them to antiproton-nucleus scatter- 
ing and to heavy ion collisions. By just imposing G-parity arguments, like in microscopic 
models 12, 14], the RMF do not describe the experimental data 0,0, 16|. This incom- 



patibility of mean-field models with respect to G-parity symmetry has been also shown 
in recent transport studies [1^, where one had to largely decrease the antinucleon-meson 
couplings by hand in order to reproduce the empirical data. 

In this work we address this issue why the conventional RMF models do not describe 
antiproton-nucleus Dirac phenomenology. To be more specific, our studies are based on the 
non-linear derivative (NLD) model jl7|] to RMF. The NLD model describes simultaneously 
the density dependence of the nuclear equation of state and the energy dependence of 
the proton-nucleus optical potential. Latter feature is missing in standard RMF models. 
Then applying G-parity transformation it is shown that the real part of the proton and 
simultaneously the real part of the antiproton optical potentials are reproduced fairly well 
in comparison with phenomenological studies. We finally make predictions for the deepness 
of the real part of the antiproton optical potential and estimate its imaginary part at low 
energies and energies relevant for the forthcoming experiments at FAIR. 

2. NLD formalism 

The NLD approach to nuclear matter is based essentially on the Lagrangian density 
of RHD [13]. It describes the interaction of nucleons through the exchange of auxiliary 
meson fields (Lorentz-scalar, a, and Lorentz-vector meson fields co^) [l8| 

^Dirac ^mes ~l~ ^int • (1) 

The Lagrangian in Eq. ([T]) consists of the free Lagrangians for the Dirac field ^ and for 
the meson fields a and u^. The isovector meson p is not considered here, for simplicity. 

In conventional RHD the interaction Lagrangian Lint contains meson fields which couple 
to the Dirac field via the corresponding Lorentz-density operators g^'^'^a and —g^^^^'^u^ 
for the scalar and vector parts, respectively. Such interactions describe rather success- 
fully the saturation properties of nuclear matter, but they miss the energy dependence 
of the mean field. A possible solution to this problem has been proposed in [ y| w here 



the momentum-dependent phenomenological form factors were introduced. In [17| this 
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idea has been generalized in a manifestly covariant way. In particular, the symmetrized 
interaction in the NLD model is given by 



r - ^ 

'~'int c\ 



9^ 
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(2) 



The interaction between the Dirac and the meson fields has a similar functional form as in 



standard RHD 13|. However, now new operators V acting on the nucleon fields appear, 
which are the non-linear functionals of partial derivatives 



V.= exp , V-= exp ' 



A 



A 



(3) 



In Eq. (|3]) denotes a dimensionless auxiliary 4-vector and A stands for the cut-off pa- 
rameter. The latter has been adjusted to the saturation properties of nuclear matter 17 . 
In the limiting case of A — t- oo the standard Walecka model is retained. 

The NLD Lagrangian £ is a functional of not only \E', and their first derivatives, but 
it depends on all higher order covariant derivatives of and For such a generalized 
functional the Euler-Lagrange equations take the form 17 
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Contrary to the standard expressions for the Euler-Lagrange equation, now infinite series 
of terms {n — )■ oo) proportional to higher order derivatives of the Dirac field {<p = 
appear. They can be evaluated by a Taylor expansion of the non-linear derivative oper- 
ators ([3]). As shown in [13], in nuclear matter an infinite series of terms can be resumed 
exactly and the following Dirac equation is obtained 



[7^(ia'^-S^)-(m-S,)]^ = 0, 
with Lorentz-vector and Lorentz-scalar self-energies defined as follows 



8 ~^ 



The Proca and Klein-Gordon equations for the meson fields can be also derived 



8,F^'' + miu'' = -g^ 
8^8^a + mla = ^g„ 



^8 8 ^ ■ 

v|/e A + '^Y^ ^ 



(5) 



(6) 



(7) 
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with the field tensor F^'^ = d^u" — 8'^u^. The meson field equations ([7]) and ([H]) show 
a similar form as in the linear Walecka model of RHD, except of the highly non-linear 
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behavior of the source terms, which generate selfconsistent couphngs between the meson- 
field equations. 

Applying the usual RMF approximation to the idealized system of infinite nuclear 
matter, the Dirac equation ([5]) maintains its original form. However, we have to distinguish 
between nucleons (A^) forming the nuclear matter and antinucleons (A^) which interact with 
the nuclear matter. For the description of antiparticles we require G-parity invariance 
of the Dirac equation and then follow the standard procedure of applying a G-parity 
transformation G = Ce*'^^^ to the negative energy states, where I2 is the operator associated 
with the 2nd component of the isospin "vector" and C is the charge conjugation operator. 
The invariance of the Dirac equation under charge conjugation requires that the auxiliary 
vector must be odd under C-parity transformation. With our choice of f'' = (1, ) for 
positive energy solutions [l7| this results in f'' = (— 1, ) for the charge conjugated Dirac 



field. This leads to the following Dirac equations for nucleons 

[7^(za^-S^)-(m-S,)]^^ = (9) 

and antinucleons 

[7^(2^'^ + S'^) - (m - S,)] = (10) 

interacting with nuclear matter, where n = and '^jq- = \E'c denote the positive energy 
and the charge conjugated Dirac fields, respectively. 

The nucleon and antinucleon self-energies entering Eqs. ([9]) and f lTOj) are the same 

ri E — m 

= S = g^uoe a , 

T,s = Qacre A . (11) 



However, note the opposite signs in the Lorentz- vector interactions in Eqs. ([9]) and f fTOj) . 
Furthermore, the single particle energies E have to be obtained from the in-medium mass- 
shell conditions which are different for nucleons (A^) and antinucleons (A^) 



E^ip) = Vp^ + + , Ej^ip) = + _ _ (12) 

The in- medium (or effective) Dirac mass in Eq. ( !T2|) is given by m* = m — E^. Note, that 
m* depends explicitly on particle momentum. Again, in the limiting case of A — 00, the 
exponential factor is equal to unity and the equations are reduced to the ones from the 
Walecka model. In the NLD model the cut-off parameter A is of natural size, i.e., of typical 
hadronic mass scale in this problem. In the following, A = 770 MeV is chosen, as in the 



original work 17 



In nuclear matter the NLD equations of motion for u and a simplify to standard 
algebraic equations 

mlu'^ = g^Pv , rnla = g^Ps (13) 



E-m 



with the corresponding density sources p^ = {^nt^ ^ n) and ps = {^^e a '^^). 
The vector density p^ is not related to the conserved nucleon density p. It has to be derived 
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Figure 1: Kinetic energy dependence of the scalar and vector Lorentz-components of the antinucleon self- 
energy in nuclear matter at densities of p = Psat 

(left), p = 2psat (middle) and p = 3psat (right) using 
the linear Walecka model (dashed lines), the linear Walecka model with rescaled couplings with the factor 
^ = 0.25 (dash-dotted) and the NLD approach (solid lines). 



from a generalized Noether-theorem 17| and reads 



9uj /7f7 



9a iTT. 



A 



A 



(14) 



The meson-nucleon couplings g^j and g^j can be taken from any linear Walecka model, 
e.g., [13|, as it has been done here. Moreover, we use the same coupling constants for both 
nucleon and antinucleon interactions. 



3. Results and Discussion 

We have applied both the NLD and the conventional linear Walecka models to nuclear 
matter at various baryon densities and also at various nucleon and antinucleon energies 
relative to matter at rest. At first, we discuss the self-energies, which are real quantities 
in RMF. Then we focus our study on the energy and density dependencies of the optical 
potential, first for in-medium proton interactions, and then for the antiproton case. 

Fig. [1] shows the Lorentz-scalar and Lorentz- vector components of the antinucleon self- 
energy in nuclear matter, and S^,, as a function of the kinetic energy at three baryon 
densities. The antinucleon kinetic energy is calculated relative to the potential depth of the 
nuclear matter at rest, i.e., E^in = Ej^ — m = a/p^ + m*^ — S^, —m. The NLD calculations 
show an explicit energy dependence for both components of the antinucleon self-energy. 
In particular, the self-energies decrease with increasing energy, for all baryon densities. 
On the other hand, with rising baryon density they increase only moderately at each 
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Figure 2: Density dependence of the scalar and vector Lorentz-components of the antinucleon self-energy 
in nuclear matter at energies of E'kin = 0.5 GeV (left), = 1 GeV (middle) and E'kin — 2 GeV (right) 
using the linear Walecka model (dashed lines), the linear Walecka model with rescaled couplings with the 
factor ^ = 0.25 [iBj (dash-dotted) and the NLD approach (solid lines). 



energy. The saturation in energy and density results from the non-hnear interaction, as 



discussed in detail in Ref. [17[. In the linear Walecka model the Lorentz- vector self-energy 
grows strongly with increasing density, while the Lorentz-scalar component saturates. Both 
components in the standard RMF are energy independent. 

For antinucleon interactions in nuclear matter the mean-field potential consists of the 
sum of scalar and vector self-energies. At vanishing momentum and at saturation density 
the linear Walecka model leads to a value of — S„ — ^ —700 MeV, which is too deep 



according to phenomenology [19|, |20|]. This feature has been always a critical problem in 
standard RMF models. Even the inclusion of non-linear self-interactions of the a field (and 
eventually of the u field) j2l| do not improve the result, since non-linear self-interactions 
become pronounced only above the saturation density. On the other hand, the NLD model 
reduces considerably the deepness of the potential at zero momentum by almost a factor 
of two. The particular difference of the potential depth at vanishing momentum between 
conventional RMF and NLD is not a trivial issue. The consequences of such an energy and 
density behavior will be discussed below when considering the optical potentials. 

As discussed in Ref. El, in order to reproduce the data from antiproton-induced reac- 
tions, the antinucleon-meson coupling constants of the Walecka model have to be rescaled 
by a factor of ^ ~ 0.2 — 0.3. Fig. [1] shows also the calculations in the linear Walecka model, 
but with rescaled couplings by a factor oi^ = 0.25 (dash-dotted curves in Fig. [1]). Indeed, 
as one can see in Fig. [T], the rescaled Walecka model [l^ reproduces the NLD results in av- 
erage. However, former results fail to reproduce the energy dependence and, in particular. 
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the density dependence, as it is demonstrated in Fig. [2l In Fig. [2] the density dependence 
(at various fixed kinetic energies) of the antinucleon self-energies is displayed. The NLD 
self-energies saturate with density and energy according to microscopic Dirac-Brueckner 
studies, as discussed in detail in \17^. In the conventional Walecka model the vector self- 
energy diverges with increasing density leading to a too strong repulsion at high densities. 
In fact, this effect of repulsive nature is softened in the rescaled model to large extend, 
however, the linear divergent behavior of the vector self-energy still remains. The NLD 
calculations agree (in average) with the rescaled Walecka model around the saturation 
density and at kinetic energies around 1 GeV only. 

The very different energy behavior of the self-energies between NLD and linear Walecka 
models influences the Schrodinger-equivalent optical potential. In general, it is extracted 
from (anti)proton-nucleus scattering and therefore it is suited for comparisons between 
theory and empirical studies. Its real part is given by 

$Ke[/opt = ±-K - + (S^ - El) , (15) 
m 2m 

where E is the energy of an (anti)nucleon with bare mass m inside nuclear matter at a 
flxed baryon density and upper (lower) sign holds for nucleons (antinucleons). At flrst we 
consider the proton-nucleus optical potential. 

Fig. E] shows the real part of the optical potential according to Eq. ( |T5|) as function of 
the nucleon kinetic energy Ekin = E — m = \/ + m*"^ + — m. The linear Walecka 
model (dashed curve) predicts the behavior of the optical potential versus energy only 
qualitatively, and strongly diverges with increasing kinetic energy of the nucleon. It does 
not reproduce the empirical saturation at higher energies. This problem is well known 
in RMF and has already attracted much attention in the past 23|. Of course, the main 
reason for such strong deviation is the missing energy dependence of the self-energy in the 
standard RMF. As discussed in detail in Ref. [13], the NLD approach resolves this issue of 
RMF models. The solid curve in Fig. |3] corresponds to the NLD calculations and describes 
the data very well. 

On the other hand, the interaction of an antinucleon at a given momentum relative to 
nuclear matter at rest is quite different with respect to the proton-nucleus interaction: the 
sign of the Lorentz- vector self-energy changes in Eq. (1151) . Therefore, in the linear Walecka 
model the real part of the optical potential is again a linear function in energy, as in the 
nucleon case, but now it diverges to — oo (see Fig. HI dashed curve). Such a prediction is in 
contradiction with calculations using dispersion relations In fact, by fltting the imagi- 
nary part of the antinucleon-nucleus optical potential to the total proton-antiproton cross 
section, its real part decreases with increasing energy. Furthermore, an existing informa- 
tion from heavy-ion collisions 0] and reactions induced by protons and antiprotons 15 
give clear evidence for a considerable reduction of the antiproton-nucleus optical potential 



with rising energy. As has been discussed in Ref. [15|, the transport theoretical description 
of antiproton-nucleus data is not possible within the conventional Walecka model, except 
if one rescales the antinucleon-meson coupling constants by a phenomenological factor of 
^ ^ 0.2. This is not compatible with G-parity arguments and suggests a strong violation 
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Figure 3: Energy dependence of the Schrodinger equivalent proton optical potential at saturation density 
Psat = 0.16 fm^^. Theoretical calculations in the linear Walecka model (dashed) and NLD approach 
(solid) are compared to Dirac phenomenology 22 1. 



of the charge conjugation symmetry in the nuclear medium [16|, which otherwise must be 
a perfect symmetry in strong interactions. 

On the contrary, the NLD calculations (solid curve in Fig. H]) predict a completely 
different behavior as compared with the Walecka model. It results in a much softer potential 
at vanishing momentum and much stronger decrease of the real part of the optical potential 
^Hcf/opt with increasing energy. 

Due to the large annihilation cross section experimental data at low energies can be 
obtained only at very low densities p ~ (0.005-i-0.02)/9sat close to the nuclear surface jl9l.l20|. 
while empirical information at saturation density is obtained by extrapolation only. At 
these low densities the NLD model leads to values of ^itUopt — —(64-50) MeV, which seem 



19, 20 



to be still too deep with respect to the data 
the NLD model predicts a rather soft potential 
data 



19l . |20| and dispersion relations 



At the density of interest p = Psat 
which is much closer to extrapolated 
(filled box at zero kinetic energy in Fig. H]). 
A comparison between our model and phenomenological antiproton-nucleus reactions at 
higher energies seems more meaningful. In fact, with increasing energy the annihilation 
cross section drops strongly and it is supposed that the antiprotons may penetrate deeper 
inside the nuclear interior, and thus densities close to psat can be tested. The second 
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Figure 4: Same as Fig. [Sj but for the antiproton case. The filled box at vanishing momentum represents 
empirical data extrapolated to saturation density [l9| . The second filled area at kinetic energies between 
1000 and 2500 MeV is taken from transport calculations on antiproton- nucleus reactions 15| . 



filled area in Fig. H] shows the empirical optical potential as extracted from the transport 



theoretical analyses in Ref. 15|, |2J]. In this energy region the comparison between NLD 



results and transport calculations (which use conventional RMF, but with largely reduced 
antinucleon-meson couplings) turns out to be fairly well. Our results are also in qualitative 
agreement with the analysis of Ref. jisf, where a strong decrease of 9let/opt with increasing 
energy is obtained. 

Interestingly, the antinucleon optical potentials DicUopt strongly differ at zero momen- 
tum between NLD and standard RMF, while in the nucleon case (see Fig. [3]) no differences 
were visible. By considering the fields at the same baryon density and at zero momentum 
one would naively expect a similar potential depth for both models. Indeed, the non- 



linear effects start to dominate above the saturation density [17|. However, the observed 
difference at zero momentum comes from the in-medium dependence of the (anti) nucleon 
single-particle energy. At fixed saturation density the energy shift, caused by the difference 
(proton-nucleus) or sum (antiproton-nucleus) of two big fields, varies strongly between the 
two models. However, small shift variations in the energy affect the NLD self-energies, due 
to their pronounced energy dependence. On the other hand, the fields of the linear Walecka 
model are not infiuenced due to their independence on energy. This feature becomes more 
pronounced with increasing density, as seen in Fig. [T] (middle and right panels) . In terms 
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Figure 5: Energy dependence of the imaginary part of the antinucleon optical potential at low density, as 
indicated. The theoretical result, extracted from the dispersion relation (see Eq. (HH)) in the NLD approach 
(solid curves) is compared to experimental data (symbols) taken from :20:]. The filled area indicates the 
model changes by varying the density from 0.005 up to 0.02 (relative to Psat)- The inserted panel shows 
the same quantity but at the fixed saturation density psat — 0.16 fm~^, again in comparison to data 
(symbols), which are extrapolated to Psat- 



of the optical potentials the interpretation is similar. In the proton-nucleus case (Fig. [3]) 
the slopes between both models at vanishing momentum are essentially the same. There- 
fore the in-medium energy shift is of minor relevance and there is no gap between both 
potentials. In the antiproton-nucleus case (Fig. H]) the gap is much more pronounced due 
to the quite different slopes between NLD and linear Walecka optical potentials. 

For a complete description of in-medium antiproton interactions also the imaginary 
part of the optical potential is needed. Since RMF does not provide the imaginary part of 
the self-energies, we estimate 3m?7opt using dispersion relation ji^ 

JmUoptip) = V / — —dp , (16) 

Jo - p^ 

where p =\p\ means the antiparticle momentum and V denotes the principal value. The 
real part 9^e[/opt is taken from the NLD model. The results are shown in Fig. [5l The insert 
in Fig. |5] shows Jmt/opt versus the kinetic in-medium energy Ekm at saturation density. 
At vanishing kinetic energy the imaginary part of the optical potential is rather large 



10 



—200 MeV) and consistent with empirical information (error bars) [20|. At high beam 
energies Jm^/opt start to decrease again, but remains rather strong. The present estimation 
seems to be in hne with the empirical study of Ref. [25j, where Jmt/opt = —135 MeV is 
essentially independent on the energy in the range from 180 MeV up to ^ 2 GeV. 

According to Refs. [lol, [2^ antiprotons penetrate the nuclear surface up to densities 
of p ^ (0.005 — 0.02)psat before annihilation. Therefore, we calculate Jmt/opt at this 
low densities, as shown in the main panel of Fig. |5l The filled area indicates the model 
calculations for matter densities p ^ (0.005 — 0.02)/9sat- The solid curve shows the model 
result at an average density of p ~ 0.01 psat- As one can see the NLD model reproduces 
the data H fairly well also at these low densities. 



4. Summary and Outlook 

In summary, the NLD model, which incorporates on a mean-field level non-linear ef- 
fects in baryon density and simultaneously in single-particle energy, has been applied to 
nucleon and antinucleon interactions in nuclear matter. We have shown that due to the 
explicit energy dependence of the self-energies the proton-nucleus optical potential is very 
well reproduced. At the same time, the NLD model predicts a much softer real part of 
the antiproton optical potential at low energies as compared to the Walecka model. We 
also find a strong decrease of the optical potential with increasing energy. These results 
are remarkably consistent with available information from reactions involving heavy ion 
and (anti)proton beams and other studies based on dispersion-theoretical approaches. A 
comparison with the conventional Walecka model has shown that the main effect respon- 
sible for a description of the in-medium (anti)nucleon optical potential originates from the 
energy dependence of the mean-field, which is absent in standard RMF models. We further 
estimated the imaginary part of the antiproton optical potential within the NLD model 
using dispersion relation. The results were in qualitative agreement with the low density 
data and empirical extrapolations at saturation density. We, therefore, conclude that RMF 
models may remain a very useful theoretical tool for the description and analysis of the 
antinucleon interactions in nuclear medium. 

The energy dependence of the real and imaginary parts of the antinucleon optical 
potential, studied in this work, is expected to be important at energies relevant for the 
PANDA experiment at FAIR. The nuclear compression due to the strong attractive antinu- 
cleon mean-field, which significantly differs between the Walecka and the NLD models, and 
also the very different energy behavior between them will affect the dynamics of antiproton- 
nucleus reactions. Thus, we expect various observable phenomena at FAIR as important 
probes for the NLD predictions. For instance, the fragmentation of the excited and radially 
expanded residual nuclei, where the energy transferred to radial expansion is expected to 
depend on the degree of compression. Strangeness production of s = —1 and especially 
s = —2 hyperons, such as cascade (S) particles, is expected to be medium dependent, in 
particular close to threshold energies. The associate formation of single-A and, in partic- 
ular, of double-strange hypernuclei is thus supposed to be also model dependent. As an 
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outlook we conclude the importance and relevance of our results for the future activities 
at FAIR. 
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